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Abstract 

We study the impact of the use of the power series expression for the primordial tensor spec- 
trum on parameter estimation from future direct detection gravitational wave experiments. The 
spectrum approximated by the power series expansion may give large deviation from the true (fidu- 
cial) value when it is normalized at CMB scale because of the large separation between CMB and 
direct detection scales. We derive the coefficients of the higher order terms of the expansion up 
to the sixth order within the framework of the slow-roll approximation and investigate how well 
the inclusion of higher order terms improves the analytic prediction of the spectrum amplitude 
by comparing with numerical results. Using the power series expression, we consider future con- 
straints on inflationary parameters expected from direct detection experiments of the inflationary 
gravitational wave background and show that the truncation of the higher order terms can lead 
to incorrect evaluation of the parameters. We present two example models; a quadratic chaotic 
inflation model and mixed inflaton and curvaton model with a quartic inflaton potential. 
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I. INTRODUCTION 



Inflation [1] is a successful paradigm not only for solving the horizon and flatness prob- 
lems, but also for explaining the origin of density perturbations in the Universe. Inflation 
predicts adiabatic and almost scale-invariant primordial fluctuations, which are in excellent 
agreement with current observations such as cosmic microwave background (CMB) and so 
on. However, no direct evidence of inflation has yet been found. During inflation, the grav- 
itational waves could also be produced whose detection can give a direct evidence of 
inflation and would be a key test of inflation. 

Early detection of the inflationary gravitational wave background may be achieved 



through its unique signature in the polarization of the CMB [3|, |4[. The ongoing satellite 
mission, Planck [5j, can detect such indirect signal of gravitational waves if the tensor-to- 
scalar ratio is r > 0.05. The next-generation experiment, such as CMBpol |6] and Cosmic 
Origins Explorer (COrE) Q], are designed to reach r ~ 10~ 3 . Moreover, the direct detection 
may be possible with space-based laser interferometers such as the DECi-hertz Interferom- 
eter Gravitational wave Observatory (DECIGO) {si, El and Big-Bang Observer (BBO) 10]. 
which would provide independent information about inflation. 

While CMB polarization experiments observe large-scale gravitational waves (k ~ 
(9(0.001) Mpc -1 ), space-based laser interferometers measure gravitational waves at milli- 
hertz frequencies {k ~ (9(10 13 ) Mpc -1 ). This millihertz frequency band is the most prospec- 
tive region for direct detection of the inflationary gravitational wave background. The 
detection becomes easier at lower frequencies, since interferometer with longer arms can 
obtain larger displacement signals by gravitational waves. On the other hand, frequencies 
below a millihertz would be contaminated by the gravitational wave background generated 



from white dwarf binaries 



We should note that there are also many other mechanisms which may generate a gravita- 
tional wave background around the millihertz frequency, such as preheating 12Hl5|. bubble 
collisions during a first-order phase transition la], self-ordering scalar fields following a 
global phase transition 20j, |21|, s econd order effects from enhanced scalar perturbations 



22 



25], topological defects [26|-|32| . supernova explosions of population III first stars 
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351 ]. gamma-ray bursts [36], and so on. However, their amplitude and frequency strongly 



depend on their unknown physics, so the millihertz band is still a window to search for the 



2 



inflationary gravitational wave background. In this paper, we focus on the gravitational 
wave background from inflation and do not consider other sources which may contaminate 
the millihertz band. 

The large difference between CMB and direct detection scales means that these two types 
of observations enable us to look at different periods of inflation, which would greatly help to 



investigate the inflaton potential 
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43| . However, we should carefully choose the method 



to connect the two separate scales. A common method is to use a power-law extrapolation 
from CMB scales to direct detection scales for describing the primordial tensor spectrum. 



Yet recent works 



44H46| have pointed out that the Taylor expansion around the CMB scale 



is no longer valid at the direct detection frequency and it causes an incorrect estimation of 
the amplitude of the inflationary gravitational wave background. 

One way to avoid the wrong estimation of the spectrum is to resort to a full numerical 
calculation to obtain the gravitational wave spectrum. However, the power-law extrapolation 
is much simpler and easier than the numerical method and, in principle, its precision can 
be improved by including higher order terms in the Taylor expansion as much as possible. 
In this paper, we derive the slow-roll expression for the primordial tensor power spectrum 
up to the sixth order in the Taylor expansion and examine how much the inclusion of the 
higher order terms improves the estimation of the amplitude at direct detection scales by 



comparing with the full numerical computation |46|. Furthermore, we discuss the impact of 
the truncation of the higher order terms in the power series expansion of the tensor spectrum 
by presenting constraints on inflationary parameters expected from future direct detection 
experiments, which is an example that such a poor estimation of the spectrum amplitude 
causes a problem. 

This paper is organized as follows: In Sec. [Til we gi ye a formula of the power series 
expression for the primordial tensor spectrum including up to the sixth order in the Tay- 
lor expansion. Next, in Sec. IIIIl we discuss whether the expression given in Sec. HJJ well 
describes the tensor power spectrum by comparing those with numerically obtained spec- 
tra. We consider two example models for the comparison, the chaotic inflation model with 
quadratic and quartic potentials. Although the quartic chaotic inflation is already excluded 
by observations such as CMB, by adding the contribution from another source of fluctua- 
tions such as the curvaton, the quartic inflation model can be allowed due to the existence 
of the curvaton fluctuations, which is sometimes called mixed inflaton and curvaton model 
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501 ] . Note that this kind of mixed scenario can give sizable tensor-to-scalar ratio as 



well as large non-Gaussianity, which might be interesting from the viewpoint of near future 
observations. In Sec. IIV} we give expected constraints on the inflationary parameters for 
the above mentioned two models. In passing, we discuss to what extent the truncation of 
the tensor spectrum expression at some (lower) order leads to incorrect evaluation of the 
inflationary parameters. Finally, we conclude in Sec. |V] 



II. SLOW-ROLL FORMALISM AND POWER SERIES EXPANSION 

In the standard picture of the early universe, a scalar field (f>, the inflaton, drives super- 
luminal cosmic expansion, the inflation. The equation of motion for <p is given by 

4> + 3H<p+V'(0) = O, (1) 

where the dot and prime denote the derivative with respect to t and 0, respectively. The 
dynamics of inflation is often characterized by the slow-roll parameters. In this paper, we 
work with the slow-roll parameters which are defined in terms of the inflaton potential V 
and its derivatives as (m] 



Vv = Ml, 
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V'V" 
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R v 2 vw 

n-3 — yt 6 

a v - M Pl v3 , 



where the subscript (n) denotes the n-th derivative with respect to (p. Inflation lasts as long 
as e\r, \i]v\ 1, which are called the slow- roll conditions, and it ends when this condition 
is violated, max{ey(0 en d), w(0end)} = 1- In the slow-roll limit, the evolution of the Hubble 



parameter H(t) is given by H 2 ~ V/(3M^), where M P1 = l/V8nG is the reduced Planck 
mass. The duration of inflation is characterized by the e-folding number, N = ln(a en d/a), 
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which can be rewritten in terms of the potential, 

N 1 ' 



v_ 



(3) 



Within the slow-roll approximation, the primordial power spectra of scalar and tensor 



perturbations are given by [52|, [53] 



V s ~ [1 - {2C + l)e H + C VH ] 



H 4 



16tt 2 M^ H> 2 



k=aH 



P T ~[l-(C + l)e H ] s 



7T 2 Mpj 



# 2 



k=aH i 



(4) 
(5) 



where C = —2 + In 2 + 7 ~ —0.73 with 7 being the Euler constant and e# and rjn are the 
Hubble slow-roll parameters, e# = 2M^ l (H l / H) 2 and 77^ = 2Mp Y H" / H . Hereafter, we only 
consider the leading order for the slow-roll parameters. Then the power spectra are given 
in terms of the inflaton potential as 

1 V 3 



Vs 



Vn 



l2n 2 M^ V' 2 



k=aH 



3n 2 M^ 



V\k=aH- 



(6) 
(7) 



They are evaluated at the moment when each Fourier mode k crosses the Hubble horizon, 
as indicated by the subscript "fc = aH." It is often assumed that the values of V and its 
derivatives evolve so slowly during inflation that the spectra can be parametrized by using 
the Taylor expansion in terms of the logarithm of the wave number, 



V T {k) =P T *exp 
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where the coefficients are the parameters characterizing a deviation from the scale-invariant 
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spectrum, 

n T (k) 

7t(/c) 
S T (k) 



>T 



(k) 



dlnV T {k) 
dink ' 
driTik) 
dink ' 
daxik) 
dink ' 

d/3 T (fc) 
(iln/c ' 

d^ T {k) 
dink ' 

din A; 

The expression for the scalar power spectrum is the same except that the coefficient of the 
first term is (ris* — 1). The subscript * denotes quantities evaluated at the pivot scale, which 
is commonly taken to be the scale of the CMB, k* = 0.002Mpc _1 . The coefficients can be 
given in terms of the slow-roll parameters as 

n T {k) ~ -2ey 

a T {k) ~ -4e v [2e v - rj v ], 

f3 T (k) ~ -4e y [16e y + 2 V 2 V - Ue v rj v + 

7T (A;) ~ -4e y [192e y - 236e^ y + 72e y ?? y _ 4^ + 2 2e y £ y - - ^ 

*r(ife) ~ -4e y [3042e y - 4810e y r/ y + 2280e y r/ y - 328e y ?? y + 8?? y 

+500e 2 v e v ~ 324e vVv & + 33t## + ~ 32eya y + 1%4 + 41 

9 T (k) ~ -4e y [61440e y - 117840e y ?7y + 75200e y r? y - 18272e y r? y + 1408ey?? y - 16r/ y 
+ 12840e y £ y - 12596e yW £ y + 3000e v vUv ~ ^vUv + 408e y £ y ~ 94 W £ y 
-948e y( x y + 648e yW a y - 77 V 2 v a 3 v - 25£ y a y + 44e y r y - 16 W r y - 

(10) 

The amplitude of the tensor perturbation at the CMB scale is often parametrized by the 
tensor-to-scalar ratio: 

= 16e y *. (11) 



As we will show in the next section, the inclusion of the higher order terms in the Taylor 
expansion up to the 6th order seems to be in very good agreement with a full numerical 
calculation, which indicates that the above expressions would be precise enough to give 
correct tensor power spectra for many inflation models and can be used for parameter 
estimation from observational data. 
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III. OVERESTIMATION OF THE TENSOR POWER SPECTRUM 



In most works, it is common to simply adopt the power-law extrapolation from CMB 
scales to direct detection scales for describing the gravitational wave background spectrum. 
However, as we will show below, such a power-law extrapolation may not be valid and lead to 
an incorrect estimation of the spectrum amplitude at direct detection scales. In Figure [U we 
show the gravitational wave spectra calculated using the Taylor expansion truncating at some 
order and the one obtained from full numerical computations. Here we consider quadratic 
(0 2 ) and quartic (0 4 ) chaotic inflation models. For the </> 4 model, we consider a mixed inflaton 



and curvaton model where fluctuations from the curvaton |54H56| also contribute to cosmic 
density perturbations. This is because the quartic chaotic inflation model predicts too large 
tensor-to-scalar ratio which is already excluded by current observational data. In addition, 
the curvaton model can generate large non-Gaussianity, thus such a mixed model would be 
interesting to investigate since it can produce both sizable gravitational wave amplitude and 
large non-Gaussianity detectable in the near future observations 1 . 

As seen from the figure, the power series expression overestimates the amplitude of the 
spectrum because of the large separation between the two scales. The spectra are plotted 
using Eq. flSJ) by truncating the Taylor expansion at each order, respectively. The exact 



spectrum, which is obtained from a numerical calculation (46], is also plotted for comparison. 
The truncation of the higher order terms in Eq. (jSJ) is the cause of the overestimation because 
the contribution of the higher order terms is non-negligible as they are boosted by the n-th 
power of ln(fco.2HzA*) — 38.7, even though the coefficients of the n-th terms are suppressed as 
e n . The overestimation of the spectrum amplitude can be avoided if the slow- roll parameters 
are much smaller than [ln(/co.2Hz/fc*)] _1 — 2.58 x 10~ 2 , but this is not the case, in particular, 
for chaotic inflation models. 

Table H] lists the values of Vt at direct detection frequency (/ = 0.2Hz) for cases of the 
truncation at each order in the Taylor expansion. The degree of overestimation compared 
to the numerical result is presented in percentage. We also list the values converted to the 
density parameter of the gravitational wave background, Oqw = (d>PGw/dlnk)/p Ci o 
where p Cj o = SM^H^ is the critical density of the Universe today and pew is the energy 



1 Mixed inflaton-curvaton models have been studied in 
details. 



47 50] and we refer the readers to these papers for 
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FIG. 1: Comparison between the exact (numerically obtained) spectrum and the spectra approxi- 
mated by truncating the Taylor expansion after the first, second, third, fourth, fifth or sixth order 
terms in Eq. ([8]). The primordial tensor spectra Vt are plotted against frequency, / = k/2iT. The 
upper panel shows the case with the quadratic chaotic inflation model, V = m 2 (f) 2 /2. The bottom 
panel shows the case for a mixed inflaton and curvaton model with the quartic inflaton potential, 
V = A(/> 4 /4. The fraction of the curvaton contribution is fixed requiring that the tensor-to-scalar 
ratio is r = 0.1. 
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density of the gravitational waves. The primordial tensor spectrum Vx{k) can be converted 
to the present-day density parameter by using the transfer function as 



1 / A; x 2 



^gw — ^) V T {k)Tl{k). (12) 



The transfer function is given by 



4/3 

Tm = (i - n A y [ ) ( ) ( ^~y 2 J (i + + ^O, («) 



g*(Thc)\ ( g*so 



9*o J \9*s(T hc )J \y/2(kr ) 
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The effective 



where r = 2if ~ 1 , x eq = k/k eq and k eq = t^ 1 = 7.1 x 10 _2 r2 m /i 2 Mpc _1 
number of degrees of freedom is given as g*(T hc ) = g* s (T hc ) = 106.75, when the contribution 
from the relativistic standard model particles are taken into account. The values at present 
are g*o = 3.36 and g* s o = 3.90. If we assume the cosmological parameters to be Q A = 0.734, 
Q m h 2 = 0.1334, h = 0.710 (taken from the WMAP 7- year mean values 59]), the amplitude 
of the primordial spectrum at the direct detection frequency / = 0.2 Hz is given by 

^gw,o.2 Hz = 1-36 x 10- 6 P T , . 2Hz . (14) 

Figure [2] plots the amplitude of the present-day tensor spectrum at direct detection scale 
(/ = 0.2Hz) in terms of Ogw for different order truncation for the above mentioned two 
models. From the figure, we see that the inclusion of the higher order terms improves the 
overestimation significantly. 

In figure [31 we show how much the overestimation of the amplitude affects determination 
of tit- One may try to determine the tilt of the spectrum if the amplitude of the 
gravitational wave background is determined by both CMB and direct detection experiments. 
However, truncation of the higher order terms would yield wrong value of riy. The values 
listed in table [I] and plotted in figure [3] are estimated with Eq. (jSJ) truncating at each order, 
with the assumption that Vt is determined at both the CMB k* and direct detection scales 

Below we give some detailed discussion for models considered here: the quadratic chaotic 
inflation and a mixed inflaton and curvaton model with a quartic inflaton potential. 



A. d> 2 model 



In the case of the chaotic inflation with a quadratic potential, 



V = \m 2 4>\ (15) 
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FIG. 2: Comparison of the overestimation for different order truncation. The vertical axis shows the 
amplitude of the gravitational wave background spectrum Ogw a t the direct detection frequency 
/ = 0.2Hz. The points represent values calculated by Eq. (JSj) truncated at each order. The straight 
lines correspond to the exact values obtained from the numerical calculation. 
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FIG. 3: Comparison of the values of tit estimated with Eq. ([8]) truncated at different order. 
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TABLE I: Summary of the amplitude of the primordial tensor spectrum Vt, the present-day 
density parameter of gravitational wave background Qqw an d percentage of overestimation due 
to the Taylor expansion, which are all evaluated at the direct detection frequency / = 0.2Hz. The 
values of re^ evaluated with the truncated expression of the spectrum are also listed. 



the slow-roll parameters are given as 

tV = " V = (16) 



and Eq. ([3]) gives 



From Eq. fl6]) we obtain 



f 2 _ 3 _ T 4 _ (X> _ n 



^=4^-2- < 17 > 



V s ^ r-^rm 2 fli 4 , (18) 

6 967r 2 M|! r ' v ; 
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and from Eq. 

p - * zmr 2 ^ (19) 

The spectrum is calculated assuming that the e-folding number corresponding to the 
CMB scale is N* = ln(a end /a^) = 60, which gives ey* = 8.26 x 10~ 3 and n T * ^ —1.65 x 10~ 2 , 
a T * ^ -2.73 x 10" 4 , /3 T * ~ -9.03 x 10" 6 , 7 T * ~ -4.48 x 10~ 7 , 5 T ^ -2.96 x lO" 8 , 
6t — —2.45 x 10~ 9 . The mass of the inflaton field m = 1.53 x 10 13 GeV is determined 
to satisfy the normalization of the scalar perturbations, Vs* = 2.43 x 10~ 9 , which gives 
Vt* — 3.21 x 10~ 10 . In this model, the tensor-to-scalar ratio and the scalar spectral index 
are r = 0.132 and n s = 0.967, respectively. 



B. <p 4 model with the curvaton 

If we consider the chaotic inflation model with a quartic potential 

1 



the slow-roll parameters are given as 



V = ^A0 4 , (20) 



Mi 



Ml 3 

-3 _ ooa -^Pl _ 3 3 

y ~ 0^ ~ 4 y ' 
r£ = $ = 0, 



and Eq. (j3J) gives 

* = (22) 
8M 2 t 2 v ; 

We again take the value of the e-folding number as AT* = 60, which gives ey* = 1.63 x 10~ 2 

and r = 0.26. This large tensor-to-scalar ratio is already excluded by current observational 

constraints, but it can be avoided by introducing the curvaton fluctuations. 

In the curvaton scenario, the fluctuations in the curvaton field a produce the scalar 



perturbations, which results in a different expression for the tensor-to-scalar ratio 47H50]. 



Since our interest is in the case where the detection of the gravitational wave background is 
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possible, we assume here that the tensor-to-scalar ratio is r = 0.1. In this case, fluctuations 
both from the inflaton and the curvaton contribute to the primordial curvature perturbation. 
In this model, the scalar power spectrum is given by 

V s = Vf + V s a) = (1 + a)vf , (23) 

where Vg and V s ^ are the contributions from the inflaton and the curvaton, respectively. 
a represents the ratio of the curvaton power spectrum to the inflaton one at the reference 
scale, i.e., a = Vg /Vg . For the 4 potential, Vg is given by 

<pf ~ 1 — r \ ( f) 6 . (24) 

The scalar spectral index is also modified as 

4e y - 2rj V 

n s = 1 - 2e v — . (25) 

1 + a 

Although the expressions for the scalar perturbation quantities are modified in this kind 
of mixed models, the formulae for the tensor perturbation spectrum Vt and the parameters 
for its scale dependence, n T , a T , (3 T , j T , <5 T , &t, are not modified from the usual inflationary 
predictions without the curvaton. The tensor spectrum in the 4 chaotic inflation model is 
obtained from Eq. (J7|) as 

T ^Jm^ 1 - (26) 



Since the scalar power spectrum is modified as in Eq. ( 123j) . the tensor-to-scalar ratio is given 
by 

r = T~~r ■ 27 
1 + a 

Assuming N* = 60, we obtain n Tic ~ -3.25 x 10~ 2 , a T * ~ -5.29 x 10" 4 , f3 T * ~ -1.72 x 
10~ 5 , 7t* ~ -8.39 x 10" 7 , 5 T ^ -5.46 x 10~ 8 and 6 T ~ -4.43 x 10^ 9 . Given the value 
ey± = 1-63 x 10~ 2 , our assumption of r = 0.1 corresponds to a* = 1.6. The normalization 
of the scalar perturbations, Vs* = 2.43 x 10~ 9 , is used to determine A = 5.94 x 10 -14 , which 
gives Vt* = 2.43 x 10~ 10 . With this setup, the spectral index for the scalar perturbation is 
n s = 0.961. 

Here we briefly comment on non-Gaussianity in this scenario. Usually non-Gaussianity of 
density fluctuations is represented by so-called non-linearity parameter /nl, which charac- 
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terizes the size of 3-point function or bispectrum 2 . Since the standard single- field inflation 
model predicts very small values of /nl as /nl "C (9(1) , if the values of /nl is found to be large 
in the future, it indicates that we need another source of density fluctuations other than the 



inflaton. As another mechanism of density fluctuations, the curvaton model 54 -[|36] has been 
intensively investigated, and in particular, this model can generate large non-Gaussianity. 
Even if fluctuations from the inflaton also contribute to the density fluctuations in the Uni- 
verse, as far as the curvaton also generates some fraction of the fluctuations, /nl &1 can be 
large. Furthermore, large tensor-to-scalar ratio is also possible in this model, which can be 
detectable in the near future. Note that, when the curvaton is the only source of density 
fluctuations, which is usually assumed in many works, the tensor-to-scalar ratio becomes 
very small. However, this kind of mixed model can give sizable /nl and r. 

In this mixed scenario where local-type non-Gaussianity is generated, /nl is given by 
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60] 



/nl= (j^) 2 AT at0n) ■ (28) 

Here / N ™ rva is the one for pure curvaton model (the curvaton is the only source of density 
fluctuation). Depending on the mass, the decay rate and the initial amplitude of the curvaton 
field, y^ rvaton ) can be very large. Thus, by tuning these parameters, the case of a* — 1.6 
(and r = 0.1), which is assumed in this section, can also give large /nl- Hence, once the 
gravitational waves and (large) non-Gaussianity are detected, this kind of scenario would be 



worth investigating carefully 61]. 



IV. IMPACT ON PARAMETER ESTIMATION 

Now, in this section, we study the influence of the poor estimation of the gravitational 
spectrum amplitude when one adopts the Taylor approximation truncated at some order. If 
direct detection determines the amplitude of the inflationary gravitational wave background, 
one may try to extract information on the inflaton potential and the e-folding number [41 1 

66|. 



by combining observations of CMB [62H64J and other complementary experiments 
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However, as shown in the previous section, when the power series expression of the spectrum 



2 Current constraints on local- equilateral- and orthogonal-types of /nl are (95 % C.L. 



< 74, -214 < < 266 and -410 < /X" < 6, respectively. 
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is adopted, one would overestimate the amplitude of the gravitational wave spectrum at the 
direct detection scale if one truncates the expression at some lower order. Here we present 
how such overestimation of the amplitude affects the determination of the inflationary pa- 
rameters by investigating future constraints. In this section, we again consider the models 
discussed in the previous section. Note that, in this paper, we do not consider effect of 



reheating which may change the shape o 



around the direct detection frequency 
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the inflationary gravitational wave background 



A. 4> 2 model 

If the quadratic chaotic inflation is the model realized in the nature, the inflationary 
gravitational wave background could be directly detected with fiGW,o.2Hz = 1-54 x 1CT 16 
(taken from the numerical result, given in Tabled]), which is obtained assuming = 60 
and the scalar perturbation being normalized as Vs* = 2.43 x 10~ 9 . With the power-law 
approximation, one can describe the amplitude of the gravitational wave background at 
direct detection scale as 

4 16 

fi G w,o.2Hz = 1-36 x 10- 6 P T ,exp[-2(38.7e y ,) - -(38.7e y ,) 2 - — (38.7e y ,) 3 

-|(38.7e^) 4 - ^(38.7e y ,) 5 - ^(38.7e^) 6 + ■•■], (29) 

where we have used Eqs. (jHD, (TTU]) . (fill) . <HM and ln(/c 02 HzA*) = 38.7. Notice that, from 
the above expression, the relation between Vt* and ey* can be provided once the value of 
^gw,o.2Hz is determined. The values of Vt* and ey* directly give information on the e-folding 
number N+ and the mass of the inflaton m via the following relations, 

m , = a^^, (31) 

which follow from Eqs. ( TT7|) and ( [19]) . 

In Figure H] we show parameter constraints expected from future CMB observations in 
the m — plane as well as the the values of m and N* which give f2Gw,o.2Hz = 1-54 x 1CP 16 
at the direct detection scale for several cases of the truncation in Taylor expansion at some 

n n 

order. Expected CMB constraints are derived from the Fisher matrix analysis [3|, |4j with 
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the instrumental sensitivity of Planck |5J and CMBpol jsj, taking into account the analysis 
of both temperature and polarization data up to the multipole I = 2000. The uncertainties 
on m and are obtained by transforming parameters from (ns,r,Vs*) into (m, N+) [41], 
with other cosmological parameters (h,tt b h 2 ,tt c h 2 ,T) = (0.710,0.1109,0.02258,0.088) [59 1 
marginalized over. 

Figure H] illustrates an important fact that the values of m and N* are estimated in- 
correctly when one determines the parameters from direct detection experiments using the 
power series expression with higher order terms being neglected. The lines in the m — N+ 
plane are plotted by Eqs. ( 150]) and (15T]) with parameters Tt* and ey* satisfying Eq. (T291 . 
truncated at each order. The fiducial values of m and N* are taken to be the same as in 
Sec. IIHI Neglect of the higher order terms leads to an underestimation of Vt* or an over- 
estimation of ev*, which results in an incorrect estimation of the values of m and N*. As 
seen from the figure, the deviation of the line from the true (fiducial) value becomes larger 
as the power series expansion is truncated at lower order. 

In particular, for the case of truncation at first or second order, the deviation is not 
negligible even if the error in measuring Qqw is taken into account. To present this clearly, 
we also plot the expected error in future direct detection experiments in Fig. El We assume 



that future experiments determine the value of ^gw with an accuracy of 701472 1 



/10 _16 \ 

a naw = 8.0 x lO" 18 — - , (32) 

V * 'GW / 

which is derived from the sensitivity of the BBO experiment (Detailed values for computing 



the noise spectrum are given in Ref. 69]). The Fabry- Perot type DECIGO has a similar 
sensitivity. Thus, the region within the error band would be similar to parameter space 
allowed by constraints from direct detection by DECIGO or BBO. Therefore, the use of 
power series expression of the spectrum may lead to incorrect parameter constraints from 
direct detection experiments, when one truncates it at lower order. However, if we includes 
up to the sixth order term, the estimate almost coincides with the true value. 

B. (j) 4 model with the cur vat on 

Next, we show an example of parameter estimation for the quartic potential in the pres- 
ence of the contribution from the curvaton fluctuations to the primordial scalar perturba- 
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FIG. 4: Parameter estimation for the (j) 2 model. The values of m and N* are inferred from direct 
detection of the inflationary gravitational wave background with Ogw,o.2Hz = 1-54 x 10 -16 . Each 
line represents the values derived assuming the gravitational wave spectrum is described by Eq. 
(|29j) . truncated at first, second, third, fourth, fifth and sixth order, respectively. The fiducial point 
is shown as a cross mark. The ellipses are the marginalized 2a constraints expected from Planck 
(solid) and CMBpol (dashed). 
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FIG. 5: The values of m and iV* inferred from the determination of ^gw,o.2Hz with the 2a experi- 
mental error of DECIGO/BBO. Each panel is for a different order truncation. 
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tions. In the same way as in Sec. IIIIBl we assume the tensor-to-scalar ratio to be r = 0.1. 
In this case, the amplitude of the gravitational wave background would be determined to be 
^gw,o.2Hz = 4.15 x 1CT 17 by direct detection experiments. The determination of ^gw,o.2Hz 
provides a relation between Vt* and ey+ via 

fiGW,o.2Hz = 1-36 x l(r 6 7Y„exp[-2(38.7ey,) - |(38.7e y ,) 2 - |(38.7e^) 3 

-^(38.7e y „) 4 - f (38.7 ev g 5 - ^(38.7e^) 6 + •■■], (33) 

where we have used Eqs. ([8]), ffTOl . ffT4"l) . (T2TT) and ln(/c .2HzA*) = 38.7. It can be converted 
to the information on N+ and A by 

= — - 1, (34) 
3ttV 

A = -^-Pt*, (35) 

which follows Eqs. (122]) and (126]). 

In Fig. |6l the values of A and N* obtained from the determination of ^gw,o.2Hz are shown 
for different order truncation of Eq. (|33l) . The lines in the A — iV* plane are plotted by 
Eqs. ( 134"1) and ( 135]) with parameters Vt* and e y * satisfying Eq. ( 133]) . For the same reason 
as described in the previous subsection, the truncation of the higher order terms gives an 
incorrect estimation of the values of A and iV*. The deviation from the true value is larger 
than the <fi 2 case, because of the larger overestimation of the spectrum as presented in Sec. 
IIHI In this case, the overestimation may come not only from the truncation of the higher 
order terms of the power series expansion in terms of In k, but also those of the slow-roll 
approximation. In our numerical calculation, the slow-roll parameter is ey ~ 4.56 x 10 -2 
when the mode corresponding to 0.2 Hz exits the horizon during inflation. This means 
the second order slow- roll correction in Vt,o.2Uz (see Eq. fl5])) can be a few percent around 
direct detection scales. Note that this cannot be improved even if we take into account the 
second order slow-roll correction as long as the spectrum is extrapolated from CMB scales, 
since the second order slow-roll correction is still small (ey* ~ 1.63 x 10~ 2 ) when the modes 
corresponding to CMB scales exit the horizon. 

In Fig. [71 the lines are shown with the experimental error of direct detection experiments, 
estimated by Eq. ( 132]) . The larger error than the 2 case is because of the smaller amplitude 
of the tensor spectrum due to the reduced normalization. Furthermore, since the CMB 
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constraints are obtained marginalizing over not only the cosmological parameters but also 
a* characterizing the contribution of the curvaton fluctuations, the uncertainty becomes 
larger compared to the case for the 2 chaotic inflation model without the curvaton. 

V. CONCLUSION 

Inflation robustly predicts a stochastic gravitational wave background with a nearly scale- 
invariant spectrum. The detection of such gravitational waves is one of the next challenges in 
observational cosmology. If both CMB polarization and direct detection experiments achieve 
the detection, the independent information from the two different scales would provide a 
breakthrough in understanding the underlying physics of inflation. 

Since the two different experiments measure gravitational waves at wavelengths sepa- 
rated by 16 orders of magnitude, the deviation from the scale-invariant spectrum, which 
is traditionally expressed by the power series expansion of In k, causes a large difference in 
amplitude of the primordial spectrum between two scales. The difference comes not only 
from the first order term of the power-law expansion, so-called the spectral index, but also 
from higher order terms, so-called runnings. We have presented that, in the case of chaotic 
inflation, the truncation of the running terms leads to the overestimation of the spectrum 
amplitude at the direct detection frequency. The overestimation is more prominent in the 
case where inflation predicts large slow-roll parameters. If we consider a single-field infla- 
tion model, large slow-roll parameters correspond to a large tensor-to-scalar ratio, in case of 
which we expect to detect the inflationary gravitational waves. Therefore, the overestima- 
tion of the tensor power spectrum should be carefully taken into consideration in case we 
achieve detection of the inflationary gravitational wave background. 

Furthermore, we have investigated how the overestimation affects the determination of 
inflationary parameters. We have considered parameter constraints obtainable from future 
direct detection experiments, assuming a specific form of the inflation potential. We have 
presented two examples: a quadratic chaotic inflation model and mixed inflation and curva- 
ton model with a quartic inflaton potential. In both cases, the use of truncated power-law 
spectrum causes an incorrect estimation of the inflationary parameters and it can be im- 
proved by adding higher order terms. For correct estimation of inflationary parameters, we 
need to take into account higher order terms, perform a numerical calculation, or develop a 
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FIG. 6: Parameter estimation for the (ft 4 model with the curvaton. The tensor-to-scalar ratio is 
assumed to be r = 0.1, which corresponds to direct detection with ^gw,o.2Hz = 4.15 x 10~ 17 . Each 
line represents the values derived assuming the gravitational wave spectrum is described by Eq. 
(f33l) . truncated at first, second, third, fourth and sixth order, respectively. The ellipses are the 
marginalized 2a constraints expected from Planck (solid) and CMBpol (dashed). 
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FIG. 7: The values of A and iV* inferred from the determination of ^gw,o.2Hz with the 2a experi- 
mental error of DECIGO/BBO. Each panel is for a different order truncation. 
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new parametrization of the spectrum to connect the two separate scales. 
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